VlO.l The Divergence Theorem 



1. Introduction; statement of the theorem. 

The divergence theorem is about closed surfaces, so let's start there. By a closed surfacx; 
S we will mean a surface consisting of one connected piece which doesn't intersect itself, and 
which completely encloses a single finite region D of space called its interior. The closed 
surface S is then said to be the boundary of D; we include S' in D. A sphere, cube, and 
torus (an inflated bicycle inner tube) are all examples of closed surfaces. On the other hand, 
these are not closed surfaces: a plane, a sphere with one point removed, a tin can whose 
cross-section looks like a figure-8 (it intersects itself), an infinite cylinder. 

A closed surface always has two sides, and it has a natural positive 
direction — the one for which n points away from the interior, i.e., points 
toward the outside. We shall always understand that the closed surface has 
been oriented this way, unless otherwise specified. 

We now generalize to 3-space the normal form of Green's theorem (Section V4). 

Definition. Let F(a;, y, z) = Mi + N j + Pk be a vector field differentiable in some 
region D. By the divergence of F we mean the scalar function div F of three variables 
defined in D by 
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The divergence theorem. Let S be a positively-oriented closed surface with interior 
D, and let F be a vector field continuously differentiable in a domain contatining D. Then 

(2) jjFdS = jjj divFdV 

We write d.V on the right side, rather than dxdydz since the triple integral is 
often calculated in other coordinate systems, particularly spherical coordinates. 
The theorem is sometimes called Gauss' theorem. 

Physically, the divergence theorem is interpreted just like the normal form for Green's 
theorem. Think of F as a three-dimensional flow fleld. Look first at the left side of (2) . The 
surface integral represents the mass transport rate across the closed surface S, with flow out 
of S considered as positive, flow into S as negative. 

Look now at the right side of (2). In what follows, we will show that the value of div F 
at {x, y, z) can be interpreted as the source rate at {x, y, z): the rate at which fluid is being 
added to the flow at this point. (Negative rate means fluid is being removed from the flow.) 
The integral on the right of (2) thus represents the source rate for D. So what the divergence 
theorem says is: 

(3) Sux across S = source rate for D ; 

i.e., the net flow outward across S is the same as the rate at which fluid is being produced 
(or added to the flow) inside S. 
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To complete the argument for (3) we still have to show that 

(3) div F = source rate at {x, y, z) . 

To see this, let Pq : {xQ,yo,Zo) be a point inside the region D where F is defined. 

(To simplify, we denote by (div F)o, {dM/dx)o, etc., the value of these functions at Pq.) 

Consider a little rectangular box, with edges Ax, Ay, Az parallel to the coordi- 
nate axes, and one corner at Pq. We take n to be always pointing outwards, as ^ 
usual; thus on top of the box n = k , but on the bottom face, n = — k . 

The flux across the top face in the n direction is approximately 

F{xo,yo,zo + ^z) ■ k AxAy = P{xo,yo,zo + ^z) AxAy, 

while the flux across the bottom face in the n direction is approximately 

F{xo,yo,zo) ■ - AxAy = -P{xo,yo,zo) AxAy . 
So the net flux across the two faces combined is approximately 

[P{xo,yo,zo + Az) - P{xo,yo,zo)] AxAy = (^^^ AxAyAz. 

Since the difi^erence quotient is approximately equal to the partial derivative, we get the 
first line below; the reasoning for the following two lines is analogous: 

fdP\ 

net Bux across top and bottom w — — AxAyAz; 

\9z Jo 

fdN\ 

net Bux across two side faces ~ -— AxAyAz; 

\ dy Jo 

(dM\ 

net Bux across front and back w — — AxAyAz; 

\ 9x Jo 

Adding up these three net fluxes, and using (3), we see that 

source rate for box = net Bux across faces of box 
fdM dN dP\ , , , 

Using this, we get the interpretation for div F we are seeking: 

„ source rate for box , ^, 

source rate at Pq = lim ; — = (div F)o . 

box->o volume ot box 

Example 1. Verify the theorem when F = xi + yj + zii and S is the sphere p = a . 

Solution. For the sphere, n = ; thus F • n = a, and / / F • rfS = Ana^ . 
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On the other side, div F = 3, JJJ 3 dV = 3 • ^""ti^; thus the two integrals are equal. □ 

Example 2. Use the divergence theorem to evaluate the flux of F = a;^ i + j + z^ii 

across the sphere p = a. 

Solution. Here div F = 2,{x'^ + y'^ + z^) = ip^ . Therefore by (2), 

11 F-dS = sJJJ pUV = s£ p'^-Anp-'dp = ; 

we did the triple integration by dividing up the sphere into thin concentric spheres, having 
volume dV = Anp^ dp. 

Example 3. Let 5*1 be that portion of the surface of the paraboloid z = 1 — .t^ — 
lying above the xy-plane, and let ^2 be the part of the xy-plane lying inside the unit circle, 
directed so the normal n points upwards. Take Y = yzi +xzi +xy'k; evaluate the flux of 
F across by using the divergence theorem to relate it to the flux across ^2. 

Solution. We see immediately that div F = 0. Therefore, if we let ^2 be the same 
surface as ^2, but oppositely oriented (so n points downwards), the surface 5*1 + 5*2 is a 
closed surface, with n pointing outwards everywhere. Hence by the divergence theorem, 

// FdS = 0= // F-dS - [ [ F-dS 

Therefore, since we have n = k on 52, 

// F-dS= // F-dS = II F kd5= // xydxdy 

= 0, 

by integrating in polar coordinates (or by symmetry). 
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